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Electromagnetic  penetration  into  a  matched  rectangular  waveguide  through  a 
longitudinal  slot  in  a  conducting  plane  is  considered.  The  problem  is  divided 
into  two  parts,  one  in  half  space  and  another  inside  the  waveguide,  by  use  of 
the  equivalence  principle, Ba— described  in  reference  [1].“  The  resultant  equa¬ 
tions  are  the  same  as  those  described  in .  reference -R]  The  only  difference 
is  that  the  excitation  vector  is  changed jfrom  an  incident  wave  inside  the  wave¬ 
guide  to  an  incident  plane  wave  in  half  Space.  The  magnetic  current  over  the 
slot  is  obtained  from  the  solution  of  the  matrix  equations.  The  electroma 
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20.  ABSTRACT  (continued  ) 


-  magnetic  penetration  into  a  matched  rectangular  waveguide  and  the  electro¬ 
magnetic  reradiation  into  half  space  through  a  slot  are  computed  from  the 
equivalent  magnetic  current  over  the  slot.  The  computed  results  show  that 
a  plane  wave  at  normal  incidence  and  a  resonated  slot  near  the  side  of  the 
waveguide  broad  wall  allow  orders  of  magnitude  greater  power  to  penetrate 
into  the  waveguide  than  when  the  slot  is  not  resonated. 
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I.  INTRODUCTION 


The  electromagnetic  penetration  through  an  aperture  in  a  conducting 
plane  is  widely  encountered  in  electromagnetic  engineering  and  is  an  im¬ 
portant  problem  in  the  engineering  analysis  of  electromagnetic  compatibility 
and  interference.  Here  we  consider  the  electromagnetic  penetration  into  a 
matched  rectangular  waveguide  through  a  longitudinal  slot  in  a  conducting 
plane.  At  resonance,  orders  of  magnitude  greater  power  can  penetrate  a  slot 
than  the  same  slot  not  resonated. 

A  general  formulation  for  a  two-region  aperture  problem  can  be  de¬ 
rived  by  use  of  the  equivalence  principle  [3].  First  the  operator  equations 
are  obtained  in  terms  of  an  unknown  equivalent  magnetic  current,  and  then 
they  are  reduced  to  matrix  equations  via  the  method  of  moments  [4].  The 
only  coupling  is  through  the  slot,  whose  characteristics  are  described  by 
two  slot  admittance  matrices,  one  for  half  space  and  another  for  the  inside 
of  the  waveguide.  These  two  admittance  matrices  depend  only  on  the  region 
considered  ana  are  independent  of  the  other  region.  The  slot  coupling 
is  expressed  as  the  sura  of  the  two  independent  slot  admittance  matrices, 
and  the  source  term  is  related  to  the  incident  magnetic  field.  The  equi¬ 
valent  magnetic  current  over  the  slot  can  be  found  from  the  matrix  equation. 

To  obtain  the  admittance  matrix  for  the  inside  of  the  matched  wave¬ 
guide  we  use  Green's  functions  to  represent  the  fields  inside  the  waveguide, 
triangle  functions  for  the  expansion  functions  of  the  equivalent  magnetic 
current,  and  also  for  the  testing  functions. 

When  the  distribution  of  the  equivalent  magnetic  current  is  obtained, 
the  other  electromagnetic  parameters  can  be  derived  by  known  methods  [5]. 

The  penetration  power  from  ualf  space  into  the  matched  waveguide,  or  the 


reradiation  power  into  half  space  from  the  slot,  is  a  measurement  which 
depends  quadratically  on  the  equivalent  magnetic  current  over  the  slot. 

It  can  be  computed  by  use  of  the  integral  of  Poynting  vector  over  the  slot. 
As  examples,  the  computed  results  for  the  magnetic  current  distribution, 
the  resonant  length  of  the  slot,  the  penetration  power  into  the  waveguide, 
the  reradiation  power  gain  pattern,  and  the  transmission  cross  section  of 
the  slot  are  presented. 

II.  FORMULATION  OF  THE  PROBLEM 


Fig.  1.  The  original  slot  coupling  problem. 

Figure  1  represents  the  problem  of  the  slot  coupling  between  two 
regions,  one  inside  the  waveguide  called  region  a  and  another  in  half 
space  called  region  b.  In  region  b  there  is  a  plane  wave  as  an  impressed 
source,  and  region  a  is  assumed  source  free.  Each  region  of  Fig.  1  is 
bounded  by  an  electric  conductor.  Region  a  is  closed  in  the  x  and  y  direc¬ 
tions,  and  open  to  infinity  in  the  r  direction.  Region  b  is  open  to  infinity 


The  problem  is  primarily  that  of  finding  the  tangential  electric  field 
over  the  slot  and  secondarily  that  of  finding  the  fields  throughout 


space. 


conixjctor 


(b) 

Fig.  2.  The  original  problem  divided  Into  two  equivalent 
problems: 

(a)  equivalent  to  region  a 

(b)  equivalent  to  region  b 


The  equivalence  principle,  described  in  reference  [1] ,  is  used 
to  divide  the  problem  into  two  equivalent  problems  as  shown  in  Fig.  2. 

In  region  a,  the  fields  are  produced  by  the  equivalent  magnetic  current 

M  =  n  X  E  (1) 

over  the  slot  region,  with  the  slot  covered  by  an  electric  conductor. 

Here  JE  is  the  electric  field  in  the  slot  of  the  original  problem  and  n 
is  the  unit  normal  pointing  outward.  In  region  b,  the  fields  are  produced 
by  the  incident  plane  wave  and  the  equivalent  magnetic  current  -M  over  the 
slot  region.  The  fact  that  the  equivalent  magnetic  current  in  region  b  is 
the  negative  of  that  in  region  a  ensures  that  the  tangential  component  of 
the  electric  field  is  continuous  across  the  slot.  The  remaining  boundary 
condition  to  be  applied  is  continuity  of  the  tangential  component  of  the 
magnetic  field  across  the  slot.  This  condition  gives  us  an  equation  from 
which  to  calculate  the  unknown  quantity  M. 

To  express  this  boundary  condition  in  an  equation  form,  we  denote 
the  electric  and  magnetic  fields  in  region  a  as 

E®  =  E^(M)  (2) 

=  H®(M)  (3) 

where  H^(M)  are  the  fields  produced  by  M  in  Fig.  2(a),  with  the  slot 

covered  by  an  electric  conductor.  Similarly,  we  denote  the  electric  and 
magnetic  fields  in  region  b  as 

E^  •=  E^(-M)  +  +  E^  (4) 

h'’  - 

where  E^(-M),  h'’(-M)  are  fields  produced  by  -M  in  Fig.  2(b),  and  E^, 


(5) 
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are  fields  from  the  incident  plane  wave,  with  the  slot  covered  by  an 
electric  conductor.  The  last  equality  in  (5)  is  a  consequence  of  the 
linearity  of  the  operator. 

To  satisfy  the  boundary  condition  that  the  tangential  component 
of  H  must  be  continuous  across  the  slot,  we  equate  tangential  components 
of  (3)  and  (5) ,  and  obtain 

-  H^(M)  =  -  H J  (6) 

over  the  slot.  The  subscript  t  denotes  the  tangential  component  over 
the  slot.  This  is  the  basic  operator  equation  for  determining  the  equi¬ 
valent  magnetic  current  M. 

We  next  reduce  (6)  to  a  matrix  equation  using  the  method  of 
moments  to  obtain  an  approximate  solution.  For  this,  we  define  a  set 
of  expansion  functions  (M  ,  n=l,2, . . . ,NT},  and  express  the  magnetic 

“Tl 

current  over  the  slot  as 


M  =  y  V  M  (7) 

—  n  — n 

n 


where  the  coefficients  are  to  be  determined.  Substituting  (7)  into  (6) 
and  using  the  linearity  of  the  operator,  we  obtain 


(8) 


over  the  slot.  For  the  slot  region,  we  define  a  set  of  testing  functions 
{W  ,  m=*l ,  2 , . . .  ,NT}  and  a  symmetric  product 


'  A,B> 


slot 


A  •  B  ds 


(9) 


slot 

We  take  the  symmetric  product  of  (8)  with  each  testing  function  and  use 
the  linearity  of  the  symmetric  product  to  obtain 


•  •  *  •  *  •  “  »  •  »  '  ^  •  •  .  *  , 
Je  -s’*. 


V  ,  H‘  (M  )'■ 
n  m  t  — n 


V  nc  ,  H  (M  ): 
n  m  X  — n 


=  -  W  ,  H  > 
m  t 


m  =  1,2 . NT 


Solution  of  this  set  of  linear  equations  determines  the  unknown  coefficients 
\  the  equivalent  magnetic  current  M  according  to  (7).  Once  M  is  known, 
the  fields  and  field-related  parameters  may  be  computed  by  standard  methods. 

The  above  solution  can  be  put  into  matrix  notation  as  follows.  De¬ 
fine  a  slot  admittance  matrix  for  region  a  as 


’■‘X'Vn- 


(y“i  .  [-<»__,, 


and  a  slot  admittance  matrix  for  region  b  as 


[Y^]  =  [-<w  ,  hNm  )>!„_ 

m  t  — n  NXN 


Define  a  source  vector 


and  a  coefficeint  vector 


Equation  (10)  is  equivalent  to  the  matrix  equation 


[Y^  +  y’^IV  = 


The  resultant  voltage  vector 


V  =  [Y^  + 


is  then  the  vector  of  coefficients  which  gives  M  according  to  (7) 


.%■  • 
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It  is  iraportant  to  note  that  computation  of  [Y  '  involves  oniv 


region  a  and  computation  of  [Y  ]  involves  only  region  b.  Hence,  we  have 
divided  the  problem  into  two  parts,  each  of  them  may  be  formulated  inde¬ 
pendently. 


III.  THE  SLOT  ADMITTANCE  MATRIX  INSIDE  THE  WAVEGUIDE 


Figure  3  shows  the  slot  configuration  and  the  coordinate  system 


Fig.  3.  The  slot  configuration  and  the  coordinate  system. 


where  a  and  b  represent  the  dimensions  of  the  cross  section  of  the 
waveguide.  The  slot,  which  is  a  narrow  rectangular  aperture  oriented 
along  the  z  direction,  is  assumed  to  be  in  the  broad  wall  of  the  rec¬ 
tangular  waveguide,  i.e.,  in  the  y*b  plane,  denotes  the  distance 

from  the  x*0  plane  to  the  middle  Line  of  the  slot,  and  is  called  the 
displacement  of  the  slot,  ix  represents  the  width  of  the  slot.  7 


8 


f:: 


denotes  the  length  of  the  slot.  To  apply  the  method  of  moments,  ■'  is 
divided  into  NS  subsections.  The  symbol  NT  denotes  the  number  of  the  triangle 
functions  which  can  be  set  up  on  the  base  of  NS  subsections. 

(17) 

NT  =  NS  -  1 

The  symbol  Iz  denotes  the  length  of  each  subsection. 


We  use  Galerkin's  method  to  evaluate  the  slot  admittance  matrix 
inside  the  waveguide,  and  select  triangle  functions  as  expansion  and 
testing  functions,  i.e.. 


W  =  M  =  U  T  (z),  p=l,2 . NT 

— p  — p  -z  p 


where  T^Cz)  is  a  triangle  function 


z  -  (p-l)Az 
Az 


ip(z)  - 


(p4-l)Az  -  z 
Az 


(p-l)Az  <  2  <  pAz 


pAz  <  z  <  (p+l)Az 


otherwise 


fc/.-.  tU-. 

>.• 


i  -  j 


Using  the  procedure  described  in  reference  [  2  ] ,  we  evaluate  the  admit¬ 
tance  matrix  inside  the  waveguide  as 


-  <W^, 


2ab  ^  3..  .2  ^  2 

n=0 


m=0  r  (Az)  r 
mn  mn 


-r  Az  -r  2Az 
/  „  mn  j.  „  nm 
-  4  e  +  e 


2  r  Az  +  I  (r  Az)^ 
mn  j  tnn 


-r  Az 


^  +  e  ]+^(4(r  Az  +  e  -1) 

mn 

-r  2Az  r  Az.2i, 
-  2e  (1-e  >  1 > 
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IV.  THE  SLOT  ADMITTANCE  MATRIX  FOR  HALF  SPACE 

A  good  solution  to  the  slot  admittance  matrix  for  half  space,  with 
the  formulation  and  the  computer  program,  is  available.  The  details  can 
be  found  in  reference  [6]. 

V.  THE  EXCITATION  VECTOR  FOR  A  PLANE  WAVE 

Suppose  the  excitation  is  a  plane  wave  incident  on  the  slot  from 
half  space,  as  shown  in  Fig.  4.  The  slot  is  along  the  z  direction,  y  is 


Fig.  4.  Plane  wave  incidence. 

normal  to  the  broad  wall  of  the  waveguide,  which  is  matched  at  both 
ends.  The  origin  of  the  coordinate  system  is  at  the  center  of  the  slot. 

9  is  an  angle  from  -z  direction  to  incidence  direction.  H  is  the 
magnetic  field  of  the  plane  wave  and  is  in  the  y-z  plane.  E  is  the 
electric  field  of  the  plane  wave  and  is  perpendicular  to  the  y-z  plane, 
it  represents  the  propagation  vector  of  magnitude  Itt/'X  and  pointing  in  the 
direction  of  propagation,  r  is  the  radius  vector  to  an  arbitrary  field 


point.  Let 


11 


where 


r  =  U  y  +  U  2 
—V  —z 


k  =  U  k  +  U  k 
-y  y  -z  z 


k  =  -  k  sin  9 

y 


k  =  k  cos  9 
z 


The  equiphase  plane  equation  is 


■  t  "jk  y  -jk  z 

-jk*r  y  z 


A  unit  plane  wave  is  given  by 


E  =  U  e 

—  —X 


~jk(-y  sin  9+  z  cos  9) 


H  =  -  (U  cos9  +  U  sin  9)e“^'^^"^  sin  9  +  z  cos  9) 
-  h  -y  -z 


where  n  is  the  intrinsic  impedance  of  the  medium  in  half  space. 

The  Impressed  field  used  in  (13)  is  the  tangential  magnetic 
field  due  to  the  incident  plane  wave  with  the  slot  covered  by  an  electric 
conductor.  In  this  case  the  y=0  plane  is  a  complete  conducting  plane  and 
image  theory  can  be  applied.  The  tangential  component  of  the  magnetic  field 
over  the  y»0  plane  when  the  slot  is  covered  by  a  conductor  is  just  twice 
what  it  is  due  to  same  plane  wave  in  free  space,  or 

hJ  -  2hJ°  (34) 

where  is  the  tangential  component  of  the  magnetic  field  over  the  slot 
due  to  the  plane  wave  in  free  space.  The  components  of  the  excitation 
vector  1  defined  by  (13)  are 


-V.N' 


V 

V*-.- 


where  W  is  the  mth  testing  function.  If  a  Galerkin  solution  is  used, 

“HI 

that  is,  if  W  =  M  =  U  T  (z),  it  then  follows 
— n  — n  — z  n 

=  -  <W  ,  H^> 
m  m  t 


=  .  f 


If  T  ^  /  N  „  IT  sin  ‘3  -jkz  cos  6  , 

1  2  T  (z)  U  *0  -  e  ds 


jJ  m  — z  — y  n 

slot 


2  sin  6 


X  +Ax/2  (iiri-l)Az  n 

or  f  -jkz  cos  6 

dx  T  (z)  e  dz 

J  ® 

X  -Ax/2  (m-l)Az 

o 


(in4-l)Az 

f  X  (z)  ®  dz 


(ni-l)Az 


(nri-1)  Az 


n  -ikmAz  cos  0  a  a  -1  A 

T  (z)  ®dz  =  ^ ^  ® 


(in-l)Az 


(k  cos  9)  Az 


^i  _  2Ax  sin  9  jkAz  cos  0  _-jkAz  cos  6  \  ^-jkniAz  cos  9 

1  =  =  =  (Z-e  -  e  /  e 

nk  cos  9Az 


NT  (35) 


When  9  =  0,  we  have 


1=0  ,  ni“l ,  2 , . . . ,  NT 

m 


When  9  =  — ,  we  have 


THE  EQUIVALENT  MAGNETIC  CURRENT  OVER  THE  SLOT 

We  can  compute  the  coefficient  vector  V  from  YT . .  YT .  and  I  ,  where 

ij  i: 

Y^j  is  the  admittance  matrix  inside  the  waveguide  described  by  (21),  (22), 
and  (23),  Y^^  is  the  admittance  matrix  for  half  space  described  by  reference 
[6],  and  I^  is  the  excitation  vector  from  the  plane  wave  described  by  (35), 

(36),  and  (37).  For  0  =  tt/6,  tt/A,  11/3,  tt/2,  2-n'/3,  3Tr/4,  and  5tt/6,  the  com¬ 
puted  results  'Vj^l  are  listed  in  Table  1  and  drawn  in  Fig.  5  and  Fig.  6. 

The  other  parameters  are  a  =  .02286,  b  =  .01016,  x^=  .01448,  w  =  .00076, 

X  =  .032,  m  =  n  =  20,  NT  =  9.  The  length  of  the  slot  is  adjusted  to  all 
o 

have  approximately  the  same  phase. 

From  Fig.  5  and  Fig.  6,  it  is  seen  that  the  magnitude  of  the  equiva¬ 
lent  magnetic  current  over  the  slot  has  a  maximum  value  when  the  angle  of 
incidence  plane  wave  is  equal  to  -n/l.  The  distribution  of  the  magnitude  of 
the  equivalent  magnetic  vurrent  over  the  slot  is  symmetric  with  respect  to 
the  center  of  the  slot  when  0  =  tt/2.  The  distribution  curves  of  the  magnitude 
of  the  equivalent  magnetic  ct^rrent  over  the  slot  for  0  =  Tr/3  and  2tt/3,  tt/4 
and  3tt/4,  tt/6  and  5Tr/6  are  also  symmetric  with  respect  to  the  center  of  the 
slot. 

Suppose  NT  =  9,  m  =  n  =  20,  a  =  .02286,  b  =  .01016,  =  .032.  -  1.0, 

w  *  .00076.  X  =  .01448,  and  0  =  Tr/2.  Here  NT  is  the  number  of  the  expansion 
o 

and  testing  function,  m  and  n  are  the  numbers  of  the  waveguide  modes  taken 
into  account,  a  is  the  length  of  the  broad  side  of  the  waveguide  cross  section, 
b  is  the  length  of  the  narrow  side,  is  the  wavelength  in  free  space,  is  the 

relative  permittivity  of  the  medium,  w  is  the  width  of  the  slot,  x^  is  the  dis¬ 
placement  of  the  slot,  and  0  is  the  incident  angle  of  the  plane  wave.  We  compute 


V 

.V/.-V. 

# 


Table  1.  The  computer  results  of  |  V 


MAGNETIC  CURRENT  DISTRIBUTION 


the  coefficient  vectors  of  the  magnetic  currents  for  various  slot  lengths, 

keeping  other  parameters  unchanged.  From  the  computed  results  we  find  that 

when  the  length  of  the  slot  reaches  a  special  value,  the  magnitudes  of  the 

components  of  the  coefficient  vector  obtain  their  maximum  values  and  the 

phases  of  the  components  of  the  coefficient  vector  are  almost  the  same.  We 

call  this  length  of  the  slot  "resonant  length,"  and  denote  it  by  L  .  For 

res 

the  present  case,  it  is  equal  to  0.47584.  Now  we  take  RAT,  „  ^  =  .95  L 

1,2,3  res 

1.05  where  RAT  represents  the  ratio  of  the  length  of  the  slot 

to  the  wavelength  in  free  space,  compute  the  coefficient  vectors  of  the 
magnetic  currents  and  obtain  the  results  listed  in  Table  2  and  drawn  in 
Fig.  7. 

Table  2.  The  computed  results  of  for  RAT^  ,  ,  =  (.95,  1,  1.05)  L _ 


1  .05  L 


36.40 

2 

;  .2429X10 

.2262 

.1974x10^ 

-31 

36.29 

.4046’<10^ 

.0786 

.3324X10^ 

-31 

36.21 

.5276x10^ 

-.0145 

.4362x10^ 

-31 

36.17 

.6042x10^ 

-.0676 

.5014x10^ 

-31 

36.16 

.6304x10^ 

-.0848 

.5236x10^ 

-31 

36.17 

.6042x10^ 

-.0676 

.5014x10^ 

-31 

36.21 

.5276x10^ 

-.0145 

.4362x10^ 

1 

-3( 

36.29 

.4046x10^ 

.0786 

.3324x10^ 

-31 

36.40 

.2429x10^ 

.2262 

.1974x10^ 

-31 

UTien  the  incident  angles  of  the  plane  wave  have  different  values 
and  the  other  parameters  are  kept  unchanged,  we  obtain  the  different 
resonant  lengths  of  the  slot.  The  computed  results  of  resonant  lengths 
for  9  =  tt/6,  t/4,  't/s,  t/2,  2tt/3,  3tt/4,  and  5Tt/6  are  listed  in  Table  3. 

Table  3.  The  resonant  lengths  for  different  incident  angles 


i  '  I  1 

''/4  :  't/3  j 

, _ i _ 

.  1 
_ 1 

2it/3 

mamm 

.47612 

From  Table  3,  we  see  that  the  resonance  lengths  are  little  changed 
for  different  incidence  angles.  The  resonant  lengths  are  the  same  for 
0  =  tt/6  and  5tt/6,  tt/A  and  3Tr/4,  and  it/3  and  211/3. 


VII.  POWER  PENETRATION  INTO  THE  MATCHED  WAVEGUIDE 


The  power  transmitted  through  the  slot  into  the  waveguide  is  a 
measurement  which  depends  quadratlcally  on  the  sources.  The  complex 
power  P^  through  the  slot  is  given  by 


P 

c 


E  X  H*  •  n  ds 


(38 


slot 

where  H*  is  the  conjugate  of  the  magnetic  field  in  the  slot  and  n  denotes 
the  unit  normal  to  the  broad  wall  of  the  waveguide.  Using  the  vector 
identity  BxC*A=AxB*£,  we  have 


P 

c 


j  j  H  ^  *  ii* 

slot 


(39 


Substituting  (1)  into  (39),  we  have 


This  involves  only  the  tangential  components  of  H,  which  in  region 
denote  by  For  M  we  use  (7)  and  obtain 


NT 


H^(M)  =  I  VX(^) 
n=l 


Substituting  (41)  and  (7)  into  (40),  we  obtain 


NT  NT  f 

p  =  -  y  y  V  V*  ' 

c  m  n  I J 

m=l  n=l  ( ■' 

slot 


M  •  H  (M  )  ds 
— n  —  — n 


From  the  definition  of  the  admittance  matrix  by  (11),  we  have 


mn 


JJ 

slot 


M  •  H  (M  )  ds 
— m  —  — n 


These  are  the  elements  of  the  admittance  matrix  of  the  slot  inside 
waveguide.  Since  M  is  a  triangle  function  which  is  real,  we  have 


-  ^[  M  •  H®*(M  ) 


mn 


J 

slot 


(M  )  ds 
— m  —  — n 


Substituting  (44)  into  (42),  we  obtain 

7  NT 

p  =  y  y  V  V  Y 

^  m-l  n-l  ”  ” 


This  can  be  written  in  matrix  form  as 


P  =  V  [Y  1*V* 
c 


a  we 

(41) 

(42) 

(43) 

the 

(44) 

(45) 

(46) 


Here  V  is  defined  by  (14)  and  the  tilde  denotes  transpose.  The  time 
average  power  transmitted  into  the  waveguide  through  the  slot  is  given  by 


From  (i7)  and  (46),  the  power  transmitted  into  the  waveguide  depends  on 
Che  magnetic  current  distribution  over  the  slot  and  the  generalized  slot 


admittance  matrix  inside  the  vc.veguide.  Table  4  and  Fig.  9  show  the 

powers  transmitted  into  the  waveguide  through  the  slot  for  various  slot 

lengths  under  condition  1,  where  NT=9,  n^n»20,  a  »  .02286,  b  =  .01016, 

=  1.0,  \  =  .032,  w  =  .00076,  x  =  .01448,  and  9  =  ■^/2.  It  can  be  seen 
r  o  o 

that  the  power  reaches  its  maximum  when  the  length  of  the  slot  is 
equal  to  .48191^,  which  we  can  call  resonant  length  2  of  the  slot  inside 
the  waveguide. 


Table  4.  The  power  penetration  versus  slot  length 

I—  Q 

I  P 

ll\ 

o 


Figure  8  shows  the  powers  transmitted  into  the  waveguide  through  the 

slot  for  various  slot  lengths  under  condition  2,  where  x  =  .0115,  x  =  .0125, 

o  o 

x  *  .0145,  x  ■  .0185,  X  *  .0225,  and  the  other  parameters  are  the  same  as 
o  o  o 

under  condition  1. 

Figure  9  shows  the  powers  transmitted  into  the  waveguide  through  the 
slot  for  various  slot  lengths  under  condition  3,  where  9  »  r/6,  tt/4,  ^73,  iT/2, 
and  the  other  parameters  are  the  same  as  under  condition  1. 

From  Fig.  8  and  Fig.  9,  we  see  that  when  0  »  't/2  and  x^  =  .02250,  i.e.. 


.09638 


.1928 


.2891 


,3855 


,4819 


.2104x10"^^  .8172x10“^° 


.1107x10“®  '  .1519x10“^  '  .2345x10"® 


,5783 


.6746 


.7710 


.8674 


.9638 


.5048x10 


”7 


.2669xl0"7 


.1995x10 


.1709x10“'^  1  .1581x10“'^ 


when  the  incident  plane  wave  is  perpe''dicular  to  the  slot,  the  slot  is 
placed  at  the  side  of  the  waveguide  broad  wall,  and  the  slot  length  equals 
resonant  length  2,  the  power  penetration  into  the  waveguide  through  the  slot 


•  • 


POWER 


Vr-N 


reaches  a  maximum. '.shen  the  slot  length  is  longer  than  the  resonant  one,  the 


power  penetration  through  the  slot  placed  somewhere  between  the 

middle  line  of  the  broad  wall  reaches  a  maximum. 

Similarly,  denotes  the  power  reradiated  into  half  space 

slot.  The  complex  power  P,  given  by 

n 


P,  =  M  E  •  H*  •  n  ds 
h 


slot 


M  •  H*  ds 


slot 


Here  H  is  produced  by  -M  in  half  space,  denoted  by 


H  =  =  -  H^(M) 


Substituting  (49)  and  (7)  into  (48),  we  obtain 


NT  NT  * 

Pu  =  “  I  I  V  V 

h  m  n 

m=l  n=l 


M  •  h'^*(M  )  ds 
— m  —  -n 


slot 


From  the  definition  of  the  admittance  matrix  by  (12),  we  have 


mn 


M  •  H  (M  )  ds 
-m  —  — n 


slot 


Since  M  is  a  real  quantity,  we  have 


mn 


J  J 

slot 


M  •  h’^*(M  )  ds 
— m  —  — n 


Substituting  (52)  into  (50),  we  obtain 

NT  NT  *  ,  * 

Pu  “  I  I  V  V 

h  m  n  mn 

m»l  n*l 


side  and  the 

from  the 


(48) 


(49) 


(50) 


(51) 


(52) 


(53) 


This  can  be  written  in  matrix  form  as 


P 


h 


V[Y^ 


]*V* 


(54) 


The  time  average  power  reradlated  into  half  space  from  the  slot  is  given 
by 


(55) 


Figure  10  shows  the  powers  reradiated  into  half  space  from  the 

slot  for  various  slot  lengths  under  condition  1.  The  power  P,  reaches 

ht 

its  maximum  when  the  slot  length  is  equal  to  .47188X^,  which  we  can  call 
resonant  length  1  of  the  slot  in  half  space. 

We  have  defined  three  resonant  lengths  of  the  slot:  denotes 

resonant  length  1  of  the  slot  and  represents  that  for  which  the  reradiated 
power  into  half  space  Pj^^  from  the  slot  reaches  its  maximum.  Lj.gg2  denotes 
resonant  length  2  of  the  slot  and  represents  that  for  which  the  power  pene¬ 
trating  into  the  waveguide  P^^  through  the  slot  reaches  its  maximum, 
denotes  the  resonant  length  of  the  slot  and  represents  that  the  total  power, 

P  =  P  +  P.  ,  reaches  its  maximum.  These  three  resonant  lengths  of  the 
o  ct  ht 

slot  are  different  but  are  close  to  each  other.  When  the  length  of  the 

slot  equals  the  resonant  length,  i.e.  L  =  L  ,  the  total  power  produced  by 

ir6S 

M  reaches  the  maximum  and  the  imaginary  part  of  the  power  is  equal  to  zero. 
Also,  the  amplitudes  of  the  components  of  the  coefficient  vector  reach  their 
maximum,  and  the  phases  of  the  components  of  the  coefficient  vector  are 
close  to  zero  (when  9  =  Tr/2) . 


where  H  is  the  riagnetic  field  from  K-  in  the  presence  of  a  complete 


conducting  plane,  and  H  is  the  component  of  the  magnetic  field  in  the 

m 

direction  of  at  v  due  to  -M  in  the  presence  of  the  conducting  plane. 
— m  — m  — 

Substituting  (7)  into  (57),  we  have 


H  Ki  =  y  V  [ [  -  M  •  h“ 

mm  ^  n  I  /  — n  — 

n  V 

slot 

This  can  be  written  in  matrix  form  as 


H  =  I™V 
m  m 


where 


I™  =  -  M  •  H®  ds 

n  — n  — 


Substituting  (19)  into  (60),  we  obtain 


I®  =  - 
n 


X  +Ax/2  (n+l)Az 

or  r 


T  (z)  H  dz 
n  z 


X  -Ax/2  (n-l)Az 
o 


where  H®  denotes  the  z-direction  component  of  the  magnetic  field  on  the 
z 

slot  from  K£  .  Since  the  x-z  plane  is  a  complete  conducting  one,  we  obtain 
— m 


H®  =  2H®°  (62) 

z  z 

Where  H®°  denotes  the  z-direction  component  of  the  magnetic  field  on  the 
z 

slot  radiated  by  K!l  in  free  space.  Substituting  (62)  into  (61),  we  have 
“tn 


(n+l)Az 


T  (z)  H®°  dz 
n  z 


(n-l)Az 


For  the  far  field,  we  take  K£  perpendicular  to  r  and  let  r^  -*■  «>.  At 


the  same  time  we  adjust  K'  so  that  it  produces  a  unit  plane  wave  in  the 

m 

vicinity  of  the  origin.  The  required  dipole  is  given  by 


•  4 


1  •  “jkr 

1  _  -j  j^e  m 

Ki  4rr  ^ 

m  m 


and  the  plane  wave  field  can  be  expressed  by 

-jk  •  r 

H  =  U  e  (of 

—  — m 

Here  U  is  a  unit  vector  in  the  direction  of  k  is  the  propagation 

— m  ~  ~!n 

vector,  and  £  is  the  radius  vector  to  an  arbitrary  field  point.  Substi¬ 
tuting  (65)  into  (63),  we  have 


■  •  »  •  »_ 


l“  =  -  2  T  (z)  e  ^  U  •  U  ds 

n  n  -m  — z 


From  (59)  and  (64),  the  far-zone  magnetic  field  is  given  by 


t  ! 


-jkr  ~  . , 

H  =  e  I™V 

m  4TTr 

m 


where  the  elements  are  given  by  (66),  and  V  is  given  by  (16). 

n 

For  symmetry  about  the  axis  of  the  slot,  we  may  calculate  the  field 
in  x-z  plane  without  loss  of  generality.  Since  the  magnetic  current  has 
only  a  z  direction  component,  there  is  no  y  direction  component  of  the 
magnetic  field  in  x-z  plane.  Therefore,  we  need  only  to  evaluate  the 
magnetic  field  on  x-z  plane  as  shown  in  Fig.  12.  From  Fig.  12,  we  have 


k  =  -  U  k  sin  <J)  +  U  k  cos 
— m  —X  — z 


U  *  U  cos  ij)  +  U  sin 
— m  -X  -z 


r  ■  U  z 
—  — z 


<i> 

(68) 

(69) 

(70) 

C«VI 


Fig.  12.  Evaluation  of  the  field  in  the  x-z  plane. 


where  k  is  the  wavenumber  of  the  excitation.  Substituting  (68),  (69) 
and  (70)  into  (66) ,  we  obtain 


(n+l)Az 


rin 


-  2Ax 


-jkz  cose 


sin  (p  dz 


(n-l)Az 


-2Ax  sin  (p  IkAz  cos  -jkAz  cosej).  -iknAz  cos  p 

2  {z-e  -e  )  e 

(k  cos  (t)  Az 


-8Ax  sin  ({)  sin^(^%-°^  *^)  ...  . 

2  -iknAz  cos  P 

e 


(k  cos  p)  Az 


n  -  1,2,...,  NT 


(71) 


From  (71),  we  have 


Substituting  (67)  into  (56),  we  obtain 


2~r  ~ 

^  R?(PJ  ^ 

n  m 


1  ji^r 


Here  I  is  given  by  (71)  and  (72),  V  is  given  by  (16),  and  Re^P.^^)  is  given 
by  (54)  and  (55). 

Figure  13  shows  the  power  gain  patterns  for  . /)  =  .094,  .472  under 

o 

the  condition  where  NT  =  9,  m  =  n  =  20,  a  =  .02286,  b  =  .1016,  \  =  0.32, 

o 

e  =  1.0,  W  =  .00076,  x  =  .01448,  and  G  =  tt/2. 
r  o 


IX.  THE  TRANSMISSION  CROSS  SECTION 


The  transmission  cross  section  of  the  slot  is  defined  as  that  area 
for  which  the  incident  wave  contains  sufficient  power  to  produce  the  radi¬ 
ation  field  H  by  omnidirectional  radiation  over  half  space.  For  unit 
m 

incident  magnetic  field,  we  have 

2Trr^n|H 
m  '  m' 


27Tr^  1h 
m  '  m' 


Substituting  (67)  into  (74),  we  obtain 


T.#  l!^!^ 


Where  I  is  given  by  (71)  and  (72),  and  V  is  given  by  (16), 


34 


Figure  14  shows  the  normalized  transmission  cross  sections  for 

“  =  “/4  and  ~/2  under  the  condition  NT  =  9,  m  =  n  =  20,  a  =  .02286, 

b  =  .01016,  '■  =  .032,  c  =  1.0,  W  =  .00076,  x  =  .01448,  and  =  .5. 

or  o  o 

X,  DISCUSSION 

This  report  describes  how  to  write  the  operator  equation  for  the 
equivalent  magnetic  current  over  the  slot  of  the  matched  waveguide  by 
using  the  equivalence  principle, and  obtain  an  approximate  solution  to  the 
operator  equation  by  use  of  the  method  of  moments.  Then  other  parameters, 
for  instance,  the  power  penetration  into  the  matched  waveguide,  the  power 
gain  pattern  of  the  reradiation  from  the  slot,  and  the  resonant  length  of 
the  slot,  can  be  computed  by  use  of  the  magnetic  current  over  the  slot. 

We  have  defined  three  resonant  lengths  of  the  slot,  i.e.  !!  ,  £  ., 

ir6S  TTBSX 

'res2'  length  of  the  slot  is  equal  to  power  penetra¬ 

tion  into  the  matched  waveguide  through  the  slot  reaches  a  maximum.  When 
the  length  of  the  slot  is  equal  to  reradiation  power  into  half 

space  from  the  slot  reaches  a  maximum.  When  the  length  of  the  slot  is  equal 
to  the  total  power  produced  by  the  slot  reaches  a  maximum,  and  at  the 

same  time  the  components  of  the  coefficient  vector  of  the  magnetic  current 
have  maximum  amplitudes  and  equal  phases. 

A  plane  wave  at  normal  incidence  causes  more  power  penetration  into 
the  waveguide  than  a  plane  wave  at  oblique  incidence. 

When  the  slot  is  near  the  side  of  the  broad  wall  of  the  waveguide, 
orders  of  magnitude  greater  power  can  penetrate  into  the  waveguide  than 
when  the  same  slot  is  near  the  center  line  of  the  broad  wall  of  the  wave¬ 
guide  . 

When  a  slotted  waveguide  is  in  a  broadband  freqiency  environment. 


such  as  an  electromagnetic  impulse,  the  slotted  waveguide  system  will 
possibly  resonant  at  some  frequencies  of  interest.  At  these  frequencies 
large  amounts  of  power  will  penetrate  into  the  slotted  waveguide. 

These  possibilities  should  be  considered  in  the  engineering 
analysis  of  such  a  system,  especially  for  perpendicular  incidence  on  a 
resonant  slot  near  the  side  wall  of  a  waveguide. 


APPENDIX 


I.  PROGRAM  DESCRIPTION 

The  computer  program  consists  of  a  main  program  and  three  sub¬ 
routines. 

The  main  program  is  used  to  organize  all  the  subroutines,  calcu¬ 
late  the  admittance  matrix  inside  the  waveguide,  and  control  the  input 
and  output  data.  Some  notations  in  the  main  program  are  described  as 
follows: 

A  and  B  denote  the  lengths  of  the  broad  and  narrow  side  of  the 
waveguide  cross  section,  respectively. 

XO  is  the  displacement  of  the  slot. 

W  is  Che  width  of  the  slot. 

REP  is  the  relative  permittivity  of  the  dielectric  inside  the 

waveguide • 

PAMO  represents  the  wavelength  in  free  space. 

M  and  N  are  indicators  of  the  wave  modes* 

NT  is  the  number  of  the  triangle  expansion  functions  for  the  slot. 

RAT  is  the  ratio  of  the  slot  length  to  the  wavelength  in  free  space. 

T3  represents  the  incidence  angle  of  the  plane  wave. 

VM  are  magnitudes  of  the  elements  of  the  coefficient  vector  of 

the  magnetic  current. 

VAM  are  values  for  phases  of  the  elements  of  the  coefficient  vector. 

POC  is  the  power  penetration  into  the  matched  waveguide  from  half 

space . 

POH  represents  the  reradiated  power  into  half  space  from  the  slot. 

POCW  is  the  total  power  produced  by  the  slot. 

GPO  represents  Che  power  gain  for  Che  reradiation  wave  from  Che  slot 


GPOl  is  the  transmission  cross  section  of  the  slot 

The  subroutine  YMAT  is  used  for  computing  the  admittance  matrix 
in  half  space  and  is  called  by  statement  CALL  YMAT  (LX,  DX,  DY,  DY,  YXX) 
in  the  main  program.  Some  notation  in  subroutine  YMAT  are  described  as 
follows: 

LX  is  the  number  of  the  subsections  of  the  slot,  equal  to  (NT+1) . 

DX  is  the  phase  shift  along  the  length  of  the  subsection  of  the 

slot  in  free  space. 

DY  is  the  phase  shift  along  the  width  of  the  slot  in  free  space. 

YXX  represent  values  of  the  elements  for  one  row  of  the  admittance 

matrix  in  half  space. 

The  subroutine  CGAUSS  solves  a  system  of  linear  algebraic  equations 
and  is  a  standard  subroutine. 

The  subroutine  DAIN  is  a  data  input  subroutine  and  is  called  by  the 
statement  CALL  DAIN  (NP,  A,  B,  XO,  W,  REP,  PAMO,  M,  N,  NT,  RATI,  EK,  T3, 


CC)  in  the  main  program 


n  n  n 


[I.  PROGRAM  LISTING 


(1)  Main  Program: 


THIS  IS  MAIN  PROGRAM 

COMPLEX  E»YA(10) , U1 » H ^ G3 »G2 » EX ( 1 0 ) » 639 , U19 » 029 
COMPLEX  UDrYXXdO)  »Y1(10»10)  »EX1  (  10  )  •  Y2  (  10  dO  ) 
COMPLEX  EI3*Y5(  lOdO)  »EI31i.P0ClJi-El2»E14d3(  lOv  10)  i.AS 
COMPLEX  EX2  ( 1 0  )  ,  GP 1 » POH ,  POC »  POC 1  ( 1 0  )  v  PGH ;l  ( 1 0  > 
DIMENSION  AM( 10) »ANN(10) »6P0f 10) kGFOI ( 10) 

CALL  DA  I N  <  NP  >  A ,  B »  XO » W » REP » PAMO ,  M » N .  NT » RAT  1 » EK  »■  T3 .  CC ) 

PI=3. 141593 

CC3=1S0/PI 

PEP=SQRT<REP) 

PAM=  PAMO/PEP 

PA=PI/A 

PB=PI/B 

PX=X0*PA 

A5=2«*PI/PAM 

E=  <  0  4  0 » 1 ) *A5*PEP/ ( 240*PI *A*B ) 

AS=A5*A3 


ai.* 


DO  1  KP=1»NP 
RAT*RAT1+EK*KP 
PL=PAMO*RAT 
PDX=0.5*W*PA 
Z1=PL/<NT+1 ) 

DO  10  1=1 »NT 
YA(I)=0. 

10  CONTINUE 

DO  11  J3=lrN 
J2=J3-1 

IF  (J3.GT.1)  GO  TO  12 

AN=U;«(U 

GO  TO  14 

1 2  AN=COS ( J2»PX ) / J2*SIN  <  J2#PDX ) /PA 

AN=8«AN«AN 
14  DO  15  J7»1»M 

J6=J7-1 

B5=J2»PA*J2»PA+J6#PB#J6*PB 

BA=B5~A5 

GZ2=BA*Z1*Z1 

G3=2«E*AN*Z1>K*3/  ( GZ2#GZ2 ) 

IF  (JZ.EO*!)  G3=0.5*G3 
G1=SQRT<ABS<BA) ) 

AA1=BS/A5 

C1=Z1#61 

CKG=AA1/C1 


BA 

3N=3IN(C:l ) 

CS=CDS(:C1) 

S1=SIN(0.5*C1 ) 

S2=S1*S1 

Ul=-SN+2.»<0.0»l))tcS2 
62=G3*16:«cS2*S2#CKG 
DO  16  I=1»NT 
II=MIN0(I»3) 

GO  TO  ( 1S» 19^20) » II 

H=2  .  *G3.t  ( 2  .  )|c  (  GZ2/3-.  -AA 1  >  f  CKGt  ( (:S--3  .  -  (  0 . 0  i.  1 )  *SN  )  :)«U  1  ) 
GO  TO  21 

H=G3*  (  G22/3+2 .  *AA1+CKG*  f  COB  (  2 -tiC  1  )  -3:i<CS+4  i-  (  0 . 0  •  1.  ) 

!  * ( 3*SN-SIN  <  2*C1 ) ) ) *U1 ) 

GO  TO  21 
A3=(I-1)*C1 

H=G2*<-SIN(A3)-(0.0i.1)*C0S(A3)  ) 

YA<I)=YA<I)+H 
CONTINUE 
GO  TO  15 

IF  <C1«GT.18)  GO  TO  25 
E2=EXF'<-C1) 

S1=E2-1 

E22=EXP<-2)kCl) 

GO  TO  26 

E2=0.0 

E22=0.0 

Sl=-1 

G2=G3)tcCKG 

G22=CABS<G2) 

GG2=AL0G ( 1 . E-35 ) -ALOG ( G22 ) 

DO  50  1=1 » NT 
II=MIN0<If3) 

GO  TO  (52»53f54> »II 

H=2*G3*  <  2*  <  GZ2/3-AA1 ) +CKG# ( E2-3 ) #S 1 ) 

GO  TO  55 

H=G3*<GZ2/3+2»AAl+CKG»(E22-3*E2+4)#Sl ) 

GO  TO  55 

Hl=-GG2/< 1-2.95) 

IF  (Cl.GT.Hl)  GO  TO  24 

H=G2*(EXP(-(I-3)*Cl/4)-EXP(-(I  +  l  )!(tCl/4)  )**4 

GO  TO  55 

H=0.0 

YA(I)=YAa)+H 

CONTINUE 

CONTINUE 

CONTINUE 

CT=C0S(T3) 

CT1=ABS<CT) 


IF  (CTl ,LE, lE-4>  GO  TO  ?6 
B  A  Z = 2  Hc  P I  /  P  A  M  *  C 1  *  Z 1 
B1=U/BA2*SIN(:T3)/BAZ*Z1/(  15)KPI) 

SN1=SIN< .5*BAZ) 

DO  90  1=1 f NT 

EX<  I  )=-Bl*SNl*SNl#(COS<  I*BAZ)-(0.0»  1  )i|cSIN<  I*BAZ)  ) 

GO  TO  91 

EXX=-U/60*Z1/PI 

no  98  1=1 » NT 

EX(I)=EXX 

BK=2*PI/PAM0 

nx=BK»zi 

nY=&K*U 

LX=NT+1 

CALL  YMATCLX.nXi-tiYrYXX) 

UD=<0.0» 1 )*Z1*W/(240*PI»PI ) 

no  80  1=1 *NT 

YXX(I)=Un*YXX<I) 

no  58  I=1»NT 

no  58  J=1»NT 

K=IABS<I-J)+1 

Y1 ( I . J)=YA<K)+YXX<K) 

Y3(I> J)=CONJG<YA(K) ) 

Y5<I»J)=C0NJG<YXX(K) ) 

no  59  1=1 »NT 

EX1<I)=EX(I) 

no  60  I=1»NT 

no  60  J=1»NT 

Y2<I»J)=Yl(l!-J) 

M1=NT 

CALL  CGAUSS ( Ml r Y2 r EXl f 0 , lE-9  f I5W ) 

P0C=(0.0,0.0) 

P0CU=<0,0»0.0) 

P0H=<0.0»0.0) 

DO  37  1=1 »NT 
P0C1<I)=<0.0»0,0) 

POH1<I)=<0,0»0.0) 

EI31=‘C0NJG<EX<I) ) 

DO  36  J=1»NT 
EI3=C0NJG<EX1 ( J) ) 

P0Cl(I)=P0Cl(I)+EI3#Y3<If J) 

POHl ( I ) =P0H1 < I ) +EI3#Y5  < I »  J  > 

P0C=P0C+EX1 ( I ) *PQC1 < I ) 

P0CU=P0CW+EX1 ( I )*EI31 
P0H=P0H+EX1 c I ) *P0H1 ( I ) 

POCT=REAL(POC) 

POHT=REAL(POH) 

T2=0.0 
DO  38  1=1*7 
T2=(I-1)*PI/12 
CT2=C0S(T2) 

CT3»ABS(CT2) 


>.■  V 


-V.^.V 


Kv.vN 

V 


IF  a:T3.LE.lE-4)  GO  TO  40 

BAZ1=2#PI/FAM*CT2»21 

SN11=SIN<0.5*BAZ1) 

B11=8#W/BA21#S1N<T2)/BAZ1*Z1 
DO  41  J*1»NT 

EX2< J)=-B11»SN11»SN11*<C0S( J*BAZl)-(O.Of 1)*SIN< J#BAZ1) ) 

GO  TO  42 
EXX1»-2*W»Z1 
DO  43  J=1»NT 
EX2< J)=EXX1 
GP1=<0.0»0.0) 

DO  4S  J=lfNT 
GP1=GP1+EX2  <  J ) *EX1 <  J) 

GP2=CABS<GP1) 

GPO ( I ) =GP2*GP2/ <  240*PAM*PAM*P0HT ) 

GPOl < I ) =GP2*GP2/ < . 2a8E+5*PAM»PAM ) 

DO  68  l^lfHT 
AM<I)»CABS<EX1<1) ) 

CC1»REAL<EX1(I) ) 

CC2«AinAG<£Xl(I) ) 

DCCaCC2/CCl 

ANNlaATANCDCC) 

ANN<I)sANNlXCC3 

URITE  <6r61)  A»Bf XOtUf PAMO»Mf NrNTf RATfREPfCCf T3 
FORMAT  <//2X» 'Aa" »F7,S»2Xf 'B=" »F7.S»2X» 'X0=' »F9.7»2X» 

!  'U*' »F9,7»2X» 'PAMO»' »F7.5/2X» r I3f IXf ' N*' ^ I3» IX f 'NT=' » 
!I2»lX»"RAT*'»F9.7flXf 'REP='»F5*3»lXf 'CC='fF9,7»lXf 'T3='»F9.7// 
GO  TO  204 
WRITE  (6 >62) 

FORMAT  <1X»'YA') 

URITE  <6»63)  (YA<I) »IalfNT> 

WRITE  (6f65) 

FORMAT  dXf'EX') 

URITE  <6f63)  (EX(I) f lalfNT) 

URITE  (6f69) 

FORMAT  (IXf'VMdlXf 'VAN') 

URITE  (6r63)  ( AM< I ) r ANN( I > r I>1 > NT ) 

FORMAT  <2Xf6E12«4) 

URITE  (6>202) 

FORMAT  < IX » ' POC ' » IX  r ' POH ' > IX  > ' POCU ' ) 

WRITE  <6f63)  POCfPOHrPOCU 
URITE  (6f64) 

FORMATdXf 'GPO') 

WRITE  (6»63)  <GP0<I)>I*1»7) 

URITE  (6 >66) 

FORMAT  (IX)- 'GPOl') 

URITE  (6>63)  (GP01(I).I»1>7) 

CONTINUE 


, 

*  *  .  •  i 


m 


(2)  Subroutine  YMAT: 


C  LISTING  OF  THE  SUiBROUTINE  YMAT 

SUBFiOUT I NE  YMAT  (  LX  »  DX  i.  DY  .  YXX ) 

COMPLEX  U » U1 » U2 » U3 » U4 , EXX » TC  < 1 1 ) »  TX ( 1 1 ) . YXX (10) 

DX2=DX*DX 

N=LX-1 

U=(0.»1.) 

U4=. 1666667*0 
YU=.5*riY 
YUD=YU*DX 
YU2=YU*YU 
YIJ3=.3333333*YU 
YU4= . 25*nX2+YU3*YU 
DO  16  1=1 » LX 
IP=I+1 

XU=(I-.5)*nX 

XU2=XU*XU 

XL=XU-DX 

XL2=XL*XL 

R1=<I-1?*DX 

R2=R1*R1 

RU1=1.-.5*R2 

U1=RU1+R1*< 1 1666667*R2)*U 
U2=<Rl-RUl*U)*YUti 
U3=-,5-.5»Rl*U 
EXX=:2.*<C0S(R1)-U*SIN<R1) ) 

R7=XL2+YU2 

R8=XU2+YU2 

R3=SaRT(R7) 

R4=SQRT<R8) 

AXU=XU*ALOG ( ( YU+R4 ) /XU ) 

AXL=XL*ALOG ( <  YU+R3 ) /ABS ( XL ) ) 

AYU=YU*ALOG( ( XU+R4 ) / ( XL+R3 ) ) 

C1=AXU-AXL+AYU 

C3=YU3*  <  XU*R4-XL*R3 )  + . 1666667* ( XU2*AXU-XL2*AXL+YU2*AYU ) 
C4=YU3* ( XU*R8-XL*R7 ) 

TC<IP)=(C1*U1+U2+C3*U3+C4*U4)*EXX 

AXU=XU*AXU 

AXL=XL*AXL 

Xl= , 5* ( YU* ( R4~R3 ) +AXU-AXL) 

X3=YU» ( ♦ 8333333E-1*  <  R8*R4-R7*R3 )  + . 125* ( XU2*R4-XL2*R3 ) > 

! +.125*(XU2»AXU-XL2*AXL) 

X4=R1*YUD*<R2+YU4) 

TX(IP)=(X1*U1+R1*U2+X3*U3+X4*U4)*EXX/DX 
16  CONTINUE 

TC(1'‘=TC(3) 

TX(  i:'=-TX(3) 

DO  20  J=2»LX 

JM=J-1 

JP=J+1 

YXX ( JM )  = , 5* ( TC  <  J ) -TX  <  JP )  +  ( J- . 5  >  #TC ( JP ) +TX  <  JM ) - 
! (J-3,5)*TC< JM) )+<TC< JP)-2.*TC< J)+TC< JM) )/DX2 
20  CONTINUE 

RETURN 


(3)  Subroutine  CGAUSS 


C 

C  THIS  IS  SUBROUTINE  4 

C 

SUBROUTINE  COAUSS < N » A » B » EPS » ISU ) 

COMPLEX  A(10»10) »B<10)»C»T 

NM1=N-1 

DO  10  K=1»NM1 

C=(0,0,0.0) 

DO  2  I=K^N 

IF  <CABS<A< IfK) ) .LE*CABS(C) )  GO  TO 
C=A(IrK) 

10=1 

2  CONTINUE 

IF  (CABS<C) .GE.EPS)  GO  TO  3 

ISW=0 

RETURN 

3  IF  <IO.EQ.K)  GO  TO  6 
DO  4  J=K»N 

T=A<K» J) 

A(Kf J)=A<I0* J) 

4  A<I0»J)=T 
T=B<K) 

B<K)=B<IO) 

B<IO)=T 

6  KP1=K+1 

C=l./C 
B<K)=B<K)*C 
DO  10  J=KF1»N 
A ( K  f  J ) =A ( K » J ) #C 
DO  20  I=KP1,N 

20  A(I*J)=A(I»J)-A(I»K)*A(K»J) 

10  B< J)=B( J)~A< J»K)*B<K) 

B<N)=B<N)/A(N»N) 

DO  40  K=1»NM1 
I=N-K 

C=<0.0»0.0) 

IP1=I+1 
DO  50  J*IPl,N 
30  C=C+A(If J)»B< J) 

40  B(I)=B<I>-C 

ISW=1 
RETURN 
END 


n  n  o 


(4)  Subroutine  DAIN: 


THIS  IS  DATA  INPUT  FOR  YAWAl 

SUBROUTINE  DA  I N  <  NP » A  f  B  »  XO  »  W » REP » P  AMO  i-  M » N » NT » RATI  f  EK  f  T3  » CC ) 
NP=3 

CC5:3#3«  141593/2 

T3=3. 141593/2 

A=. 02286 

B=, 01016 

X0-.0144S 

U=. 00076 

REP=1.0 

PAM0=.032 

RAT1=.4247 

EKs.0236 

M=20 

N-20 

NT=:9 

RETURN 


III.  PROGRAM  USAGE 


We  write  a  subroutine  to  input  the  necessary  data.  Keeping  the 
main  program  unchanged  and  using  different  data  input  subroutines,  we 
can  compute  different  parameters  of  the  slot.  The  example  is  as  follows: 


C 

C 

c 


THIS  IS  DATA  INPUT  FOR  YAUAl 

SUBROUTINE  DAIN  <  NP . A . B » XO  f  U  *  REP  >  FAMO » M  ?  N , NT • RATI >  EK  >  T3  r  CC ) 
NP=3 

CC=3*3. 141593/2 

T3=3. 141593/2 

A=. 02286 

B=. 01016 

X0=. 01448 

U=. 00076 

REP=1.0 

PAf10=.032 

RAT1=. 42826 

EK=, 02379 

M=20 

N*20 

NT=9 

RETURN 

END 


A=0. 02286 
M=  20  N=  20 


B=0. 01016  X0=0, 0144800  W=0. 0007600  PAM0=0. 03200 

NT=  9  RAT=0, 4520500  REP=1,000  CC=4. 7123890  13=1,5707960 


OM  VAN 

0.2075E+02 
0.50S9E+02 
0.4432E+02 
POC  POH  POCU 
0. 1573E-06 

GPO 

O.OOOOE+00 
0. 1628E+01 
GPQl 

O.OOOOE+00 

0.2036E-07 


0.3640E+02 

0.3617E+02 

0.3621E+02 

0.6027E-06 

0.7351E-01 

0.9194E-09 


0,3419E+02 

0.5272E+02 

0.3419E+02 

0.1501E-05 

0 . 2970E+00 

0.3714E-08 


0,3629E+02 

0,3616E+02 

0,3629E+02 

0.6120E-06 

0.6609E+00 

0,8265£-08 


0,4432E+02 

0.5059E+02 

0.2075E+02 

0. 1658E-05 

O.llOlE+01 

0,1377E-07 


0.3621E+02 

0.3617E+02 

0.3640E+0:. 

0.  1215E-0!: 

0.1478E+0 ; 

0.1348E-07 


<1 


wm 


i 


A=0. 02286  6=0.01016  X0=0. 0144800  U=0. 0007600  PAi'10  =  0 . 03200 

M=  20  N=  20  NT=  9  RAT=0 .4758400  REP=1.0b0  CC=4. 7123890  13=1.5707960 


OM  VAN 

0.2429E+02 
0.6042E+02 
0.5276E+02 
POC  POH  PQCU 
0.2409E-06 

GPO 

O.OOOOE+00 
0. 1641E+01 
GPOl 

O.OOOOE+OO 

0.3ia7E-07 


0.2262E+00 

-0.6759E-01 

-0.144aE-01 

0. 1708E-06 

0.7112E-01 

0.13.81E-0a 


0.4046E+02 

0.6304E+02 

0.4046E+02 

0.2331E-05 

0.2897E+00 

O.S62aE-08 


0.7S59E'01 

-0.8482E-01 

0.7856E-01 

-0. 1704E-06 

0.6520E+00 

0.1266E--07 


Q.5276E+02 

0.6042E+02 

0.2429E+02 

0.2572E-05 

0. 1098E+01 

0.2133E-07 


-0.1445E~01 
-0.6760E-0  !. 
0.2262E1-00 

0.2468E-:)9 

0. 1485E+01 

0.2S85E-07 


A=0. 02286  6=0.01016  XOsO . 0144800  W=0.0007600  PAM0=0. 03200 

M=  20  N=  20  NT=  9  RAT=0 . 4996300  REP=1.000  CC=4. 7123890  13=1.5707960 


VM  VAN 

0.1974E+02 
0.5014E+02 
0.4362E+02 
POC  POH  POCU 
0.1770E-06 

GPO 

O.OOOOE+OO 
0. 1654E+01 
GPOl 

O.OOOOE+OO 

0,2395E-07 


-0.3056E+02 

-0.3092E+02 

-0.3086E+02 

-0.3337E-06 

0.6870E-01 

0.9945E-09 


0.3324E+02 

0.5236E+02 

0.3324E+02 

0.1737E-05 

0.2824E+00 

0.4088E-08 


-0.3074E+02 

-0.3095E+02 

-0.3074E+02 

-0.8038E-06 

0.6429E+00 

0.9307E~08 


0.4362E+02 
0.5014E+02 
0. 1974E+02 

0. 1914E-05 

0.1095E+01 

0.1585E-07 


-0.3086E+02 

-0.3092E+02 

-0,3056E+02 

-0.1143E-05 

0.1493E  +  0.!. 

0.2161E-07 
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